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Abstract 



We study fluids of hard rods in the vicinity of hard spherical and cylindri- 
cal surfaces at densities below the isotropic-nematic transition. The Onsager 
second virial approximation is applied, which is known to yield exact results 
for the bulk properties in the limit of infinitely thin rods. This approach re- 
quires the computation of the one-particle distribution function and of the 
Mayer function which is greatly facilitated by an appropriate expansion in 
terms of spherical harmonics. We determine density and orientational profiles 
as well as the surface tension 7 as function of the surface curvature radius R. 
Already in the low-density limit of non-interacting rods 'y(R) turns out to be 
non-analytic at 1/R = 0, which prohibits the application of the commonly 
used Helfrich expansion. The interparticle interaction modifies the behavior 
of j(R) as compared to the low-density limit quantitatively and qualitatively. 
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I. INTRODUCTION 



A fluid of hard rods can be considered as the simplest model for nematic liquid crys- 
tals consisting of elongated molecules. In a seminal paper in 1949 Onsager showed [0] that 
the steric hard-body interactions alone can bring about an isotropic-nematic transition. Al- 
though the steric interactions already capture many of the essential features of liquid crystals 
their actual behavior is complicated by the presence of dispersion forces, flexibility, dipole 
moments, etc. But for certain colloidal systems of rod-like particles of synthetic or biological 
origin dissolved in a suitable solvent the hard rod model provides a quantitatively reliable 
effective description 0. Among them the ones, which are studied in most detail, are the 
tabac mosaic virus and the fd- virus with length (L) to diameter (D) ratios L/D of about 17 
and 150, respectively PJ. From a theoretical point of view the limit of infinitely thin hard 
rods is especially interesting because it represents one of the very few cases for which the 
exact density functional is known |T|. 

Even more than in simple liquids, which are composed of spherically symmetric particles, 
surface effects are of great importance for liquid crystals. In the absence of external fields 
the orientation of the bulk fluid is determined by its interaction with the container walls; 
this phenomenon is called anchoring ||. For the simplest case of a hard rod fluid near a 
planar hard wall, as studied theoretically by Holyst and Poniewierski , the wall induces 
parallel alignment of the nematic director. An isotropic-nematic interface also aligns the 
nematic phase parallel to the interface for large aspect ratios L/D while a non-trivial tilt 
angle arises for lower aspect ratios 0. In view of the substantial technical difficulties which 
are associated with the theoretical description on a truly microscopic scale, especially for 
curved surfaces, it is a natural first step to analyze the interface between the isotropic phase 
and a hard wall. The corresponding density and orientational order profiles near a planar 
wall have been determined by Poniewierski || in the framework of the Onsager theory, 
who also found indications for spontaneous biaxial order at the surface already below the 
bulk transition to the nematic phase. Mao et al. have compared this theory with computer 
simulations for finite aspect ratios L/D |J and have calculated the depletion force between 
planar walls or large spheres immersed in a solution of rods [|T^,|rj], using the Derjaguin 
approximation. 

In the present work we focus on the orientational and positional order as well as the 
surface tension near curved hard walls, taking into account the steric interactions between 
the rods. The curvature has an appreciable effect on the structure and the thermodynamics 
of the fluid if the radius of curvature R is of the order of the particle length L. Accordingly 
as possible applications one can think of the following systems: (i) rod-like particles confined 
to the interior of small pores within porous materials; (ii) colloidal suspensions of rods that 
contain a second, diluted, component of larger, e.g., spherical particles; (iii) membranes, 
especially vesicles, immersed in colloidal rod solutions, resembling, e.g., solutions of viruses. 
The curvature dependent surface tension and the depletion forces in case (ii) have been 
determined by Auvray and Yaman et al. [0,0 for fluids of non-interacting rods corre- 
sponding to the limit of infinite dilution. Their most surprising result was that the surface 
free energy does not contain a term linear in the curvature 1/R and that the quadratic term 
has different amplitudes for different signs of the curvature. This non-analyticity prohibits 
the application of the common Helfrich expansion. Thus for case (iii) above the effect of the 
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rods on the elastic properties of the membranes cannot be described by a renormalization 
of the bending rigidities as it is possible for membranes exposed to a suspension of spherical 
colloidal particles [0 or polymers ||15|| . In order to be able to assess the range of validity 
of the results obtained in the ideal limit of non-interacting rods, in the present paper we 
tackle the full problem including the inter-particle interactions by employing the Onsager 
density-functional theory (Sec. 0), which yields density profiles (Sec. HD and the surface 
tension (Sec. [TV] ). Our main results are summarized in Sec. [V] while technical details are 
presented in Appendices A and B. 



II. MODEL AND DENSITY-FUNCTIONAL THEORY 

Based on density-functional theory we study a fluid of hard sphero cylinders of length L 
and diameter D in the vicinity of a hard spherical or cylindrical surface of radius R (see 
Fig. |I|). In order to keep the numerical difficulties tractable we restrict ourselves to the limit 
D/L — > with R/L fixed. The number density of the centers of mass of these thin rods 
at a point r with orientation u' = (9',(f>') is denoted by p(r,u'). The corresponding grand 
canonical functional is given by 

PSl[{p(r, uj')}\ = J d 3 rdto'p{r, u') [in (47rA 3 p(r, uo')) -1-/3/1 + /3V{r, to')] 

+ PF ex [{p(r,u')}}. (1) 

Here = 1/T is the inverse temperature, p the chemical potential, V(r, u/) the external 
potential exerted by the hard wall and A the thermal de Broglie wavelength. Within the 
Onsager second virial approximation the excess free energy F ex is 

pF ex [{p(r,uj')}] = -~ J d 3 r 1 du[d 3 r 2 du' 2 p(r 1 ,uj' 1 )p(r2,uj' 2 )f(r 1 2,uj' 1 ,uj' 2 ) (2) 

with the interparticle vector r 12 = r 2 — i"i and the Mayer function /(r 12 , u)[, uj 2 ) which for hard 
particles equals —1 if the two particles overlap and zero otherwise. Onsager demonstrated 
that this approximation becomes exact in the limit D/L — > for the bulk properties M, and 



this is expected to hold also for surface quantities igjllj- For a planar surface this approach 
has turned out to be quantitatively reliable for D/L < 0.1 [Oj. 

As mentioned in the introduction, the corresponding problem with non-interacting par- 
ticles, i.e., in the presence of the hard wall interaction but with F ex = 0, has been analyzed 
by Yaman et al. |I3| , |14|| . Taking into account the inter-particle interaction increases the 
complexity considerably due to the ensuing non-locality and the high dimensional integra- 
tion in Eq. (fj). In the remaining part of this section we evaluate the expression for F ex 
by exploiting the symmetries of the density profile p(r, uo) using appropriate expansions in 
terms of spherical harmonics. In Sec. [TTT] the thermodynamically stable equilibrium profile 
is obtained by minimization of the density functional, which amounts to solving an integral 
equation. The value of the functional at this minimum is the grand canonical potential of 
the inhomogeneous fluid from which the surface tension at the curved walls is determined 
in Sec. [TV| 

The particle orientation is conveniently described in a local reference frame whose polar 
axis is that wall normal which runs through the rod center and whose y axis, in the case of 
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a cylindrical wall, is aligned with the cylinder axis (see Fig. [I]). Since we are interested in 
bulk densities pb below the isotropic- nematic transition, i.e., pb < 4.2 (DL 2 )^ 1 Q, we may 
assume that the number density p{r) = f duj'p(r, uj') as well as the orientational distribution 
measured in the local reference system (denoted by uj = (9,<f>)) depend only on the radial 
coordinate r in a spherical or cylindrical coordinate system, which allows us to make the 
following ansatz 

oo 

p S (r,uj) = ^Y. a ^ P + cos ^ ( 3 ) 
1=0 

and 

oo I 

p C (r,Uj) = Pb^Yl Plm( r ) Y lm(uj) (4) 



1=0 m=-l 



for a sphere (S) and a cylinder (C), respectively. The functions Pi and Y im are Legendre poly- 
nomials and spherical harmonics, respectively. For large distances from the wall the fluid is 
isotropic (p(r, uj) = p b /iir), so that a>i(r — > oo) = |^ i0 and pi m (r — > oo) = (47r) -1 / 2 5i m) oo- At 
a sphere the density does not depend on the azimuthal angle 0. At a cylinder the symmetries 

p(r, 9, 0) = p(r, 7T - 9, 0) = p(r, 9, -0) = p(r, 9,n-<f>) (5) 

imply that pi m = if I or m is odd and pi m = pim = p\ m for I and m even (here and in the 
following m = — m). 

The coordinates uj can be expressed in terms of the coordinates uj' corresponding to a 
frame fixed in space by a (position dependent) rotation. Therefore the angular integrations 
in Eqs. ([Xj) and (§) can be taken also over uj. However, the Mayer function is naturally 
expressed within a third coordinate system uj given by the interparticle vector r 12 . In order 
to perform the angular integrations in the following we will determine the transformation 
from uj to uj. The definition of the different reference frames is illustrated in Fig. [l]. 



For any uniaxial molecule the Mayer function can be expanded as [16] 



f(r 12 ,uJi,uJ 2 ) = ^2 fhhm(ri2)Y hm (wi)Y hm (w2)- (6) 

h,l2,m 

The solid angles UJi refer to a particle fixed reference system (see Fig. [I]) with its z axis parallel 
to the interparticle vector r 12 (and arbitrary x axis). The determination of the expansion 
coefficients 



fhi 2 m{riz) = j duj 1 duj 2 f(r 12 ,uj 12 = 0, u Xl uj 2 )Y l l m {oj 1 )Y* m {uj 2 ) (7) 

is discussed in Appendix [A|. Spherical harmonics in different reference systems are related 
via the rotation matrices D l nm [|T6|: 



Yl m (uJi) = Dlimi^h Vi, Xi)Yln{uJi), 1 = 1,2, (8) 
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where the Euler angles ^(ri, r 2 ), ^i(ri, r 2 ), and x«( r i> r 2) describe the rotation of the particle 
based axes (a)j) onto the surface normal based axes (u!i) (for the definition of the Euler angles 
see, e.g., Fig. A. 6 in Ref. [[[§]). Inserting Eqs. (|3|), (||), and (|8|) into Eq. (0) yields for the 
spherical case 

P F Z = -fH Pi + l)(2/ 2 + I)]" 172 (9) 

h,h,m 

x y d 3 ri d 3 r 2 a h (ri)a^r 2 )/ W2m (r 12 )I^*(V>i, t&, Xi^loOa, Vz, Xa)- 

One still has the freedom to fix the orientation of the ?/ axes in the different reference systems. 
If one chooses them to be all parallel to each other and perpendicular to the plane spanned 
by ri and r 2 the transformations described by the Euler angles become simple rotations 
around the y axis so that tpi — Xi — 0. With |L6[ 



^o(o,^o) = ^^iu^o) (10) 

one finds for a system of radial size C, outside of a spherical cavity of radius R 

-R+C r R+C 



PF e s x = -(47Tp 6 ) 2 i( 2 h + l)(2i a + I)]" 1 / dr x r\ \ dr 
x / dcosj f h i 2m (ri 2 )Y hm (r] 1 ,0)Yi m (T]2,0). 



2r 2 ai 1 {r 1 )ai 2 {r 2 ) (11) 



i 



The angles 7, rji, and rj 2 are those between the vectors r 1 and r 2 , r x and r 12 , and r 2 and r 12 , 
respectively, and r 12 = (r 2 + r\ — 2r x r 2 COS7) 1 / 2 . If we use r 12 instead of cos 7 as integration 
variable we finally obtain 



Y t-R+C r-R+C 

P F ex= oPbY] / dr x r x \ dr 2 r 2 a lx {r^a l2 {r 2 )w lxl% {r x ,r 2 ) (12) 
1 hM Jr Jr 



with 



w h>h {r u r 2 ) = - V / rfr 12 r 12 / ZlZ2m (r 12 )Y, im (^,0)Y, 2m (^ 2 ,0) (13) 

(2<l + lj^i 2 + 1J ^ J|n-r 2 | 



and 



r2 v< 2 yi 2 ™2 2 I > 1 2 

2 ' 1 '12 ' 2 ' 1 1 1 12 / -1 ,1 \ 

cos Vl ^^— ^ i 2 -, cos77 2 = ^— ^ ^. (14) 

2riri 2 2r 2 r i2 



An equivalent expression for a planar wall with surface area A can be derived along the 
same lines. In this case the Euler angles are the same for both particles because the direction 
of the surface normal is the same everywhere. Again one can choose tpi = Xi — and finds 

P F L/ A = lpbY] I dz x I dz 2 a h {z 1 )a h {z 2 )wf lh (z 1 - z 2 ). (15) 
1 h,h Jo Jo 
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where the interaction kernel wf il2 now depends only on one variable: 



p , s 87T 



(2Zi + l)(2Z 2 + l) 



/■oo 

/ dr 12 r 1 2fi 1 i 2m {ri2)Y hm (r],0)Yi 2 rn(r],0). (16) 



with?] = arccos(,2i 2 /ri 2 ). It can be shown that wi u i 2 (R+zi, R+Z2) = kirwf, {z%— Z2)+0(l/R) 
for zi,Z2 <C R. 

The cylindrical case is considerably complicated by the lower symmetry of p(r,u). Using 
Eqs. (0), ([4|), and (|8|) and performing the integrations over u>i, u 2 , as well as Z\ and 0i, where 
rj = (rj, Zj, 4>i), i = 1, 2, in cylindrical coordinates, one obtains 



l^limi ( r l)V>l 2 m 2 ( r 2)w/ imi / 2m2 (ri, r 2 ) (17) 



with 

/OO /*2"7r 
dz 12 / d0 12 /i 1 i 2m (r 12 )D^ i (V'i,^i,Xi)^m 2 (V'2,^2,X2)- 
Jo 

(18) 

Here if is the macroscopic height of the cylinder and Z\ 2 = z 2 — Zi, i2 = <f> 2 — 4>i, r i2 = 
(r\ + r 2 + z\ 2 — 2r x r 2 cos0i 2 ) 1 / 2 . In a rather lengthy calculation the dependence of the Euler 
angles on r 1; r 2 , Z\ 2) and i2 can be worked out by decomposing the rotations that connect 
the different reference systems into three successive simple rotations around (intermediate) 
coordinate axes (see App. A. 2 in Ref. jnj). It is helpful to use r i2 and u — r\ 2 — zf 2 as the 
integration variables which leads to 



/•OO 

(n,r 2 ) = -16n^2 / dr 12 r 12 f hhm (r 12 ) 

„ ■/In— ral 



ri— rz| 

x f" dti[(ri a -tt)((r 1 + r 2 ) 2 -«)(tt-(r 1 -^) a )]- 1/a (19) 

^ (ri-r 2 ) 2 



x cos [m(^i - ip 2 ) + rrnxi + ^2X2] d l ^ mi (r}i)d^ m2 (rj 2 ). 
Here the rotation matrices have been written as 



D l mn ^^x)=e- l ^d l m Me- inx (20) 



where the functions d l mn can be calculated by using Eq. (A. 65) in Ref. [HJ. The advantage of 
Eq. (|19| ) is that the inner integral can be evaluated without the time-consuming calculation 
of fhizm- I n these variables the Euler angles are 

r\ — r\ — u r\ — r? + u 
cos? 7i = — n cos7y 2 = — (21) 

tan^i = 2?"! A cos 771 tan-^ 2 = 2r 2 Acos?7 2 (22) 
tanxi = -2riA tan%2 = -2r 2 A (23) 
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with 



((ri + r 2 ) 2 — u)(u— (ri - r 2 ) 2 ) 



1/2 



(24) 



For the inside of a sphere (cylinder) the integration range for the radial integrals in 
Eq. (HD (Eq. (0)) has to be replaced by [0,R]. 

III. DENSITY PROFILES 



The equilibrium density profile minimizes the grand-canonical functional, i.e., it is a 
solution of jjjyuj = under the boundary condition p(r,u) — > Pfe/47r for r — >■ oo. By using 
the relation 



<H(r') _ 2/ + 1 
Sp(r,u) 2p b 



5(r — r') Pi (cos 9) 



(25) 



and exploiting the symmetry property ^ 1 z 2 ( r i, r 2) = ^/ 2 «i( r 2, n.) one finds for the spherical 
wall the Euler Lagrange equation 



An\ 3 p(r, 9) = exp 



^-/^M)-£;£ 



2Zi + 1 



P h (cos9) / dr'r'a h (r')w h i 2 (r,r') 



(26) 



With p (r,9) = (47rA 3 ) 1 exp[/3/^ — V(r, 0)] as the corresponding profile for non-interacting 
rods at the same chemical potential one finds that 



p(r,9) :=p(r,6)/p (r,6) =: £ A(r)P,(cosi 



(27) 



satisfies 



with 



p(r, 9) = exp 



p b 21 + 1 
47rr 2 



£ (cos 9)p l (r) 



^ / dr' r'a v (r')w w (r,r'). 

v J 



(28) 



(29) 



(Strictly speaking p cannot be defined by Eq. (|27j) for the forbidden orientations, for which 



both p and p vanish. Instead we define it by Eqs. ( P8|) and fl29|) in this region.) The function 
po(r,9) equals po/47r for orientations that are allowed by the hard wall and zero otherwise. 
The density po corresponding to the chemical potential p follows from the bulk limit of the 
density functional. For an isotropic fluid in a volume V one has 



tt'W = '^'o = pi, ( In X 3 p b - 1 - /3p + ^p 2 b v 



(30) 
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with 



vo 



(4tt) 



d 3 r 12 du} 1 du} 2 /(ri2, Ui, U) 2 ) 



71 



DL 



(31) 



The same equation without the last term holds for the ideal gas limit. Minimization yields 
Po = Pb exp(p(,u ). The allowed values of 6 for given r and i? are determined in Appendix [FJ. 
Thus based on the known function po{r,x = cos 8) the coefficients ai in Eq. ( p9|) can be 
expressed in terms of the coefficients introduced in Eq. (|27j): 



21 + 1 2tt 
2 ^ 



// */ — 1 



dxPi(x)P v (x)p {r,x) 



(32) 



where the integration over x can be carried out analytically for given I and /'. This allows 
one to calculate the coefficients by solving iteratively the following system of equations 
together with Eqs. fl29|) and (|32"|): 



A(r) 



2/ + 1 



dxp (x) exp 



^P;/(x)p/.(r) 



(33) 



The advantage of first seeking the solution for p instead of p is that the former function is 
smoother near the transition from allowed to forbidden orientations and hence can be better 
approximated with a limited number of Legendre polynomials. 

The cylindrical case can be treated completely analogously. With the expansion p(r, ui) = 
J2im u im(r)Yi m (uj) one obtains 



vi m (r) = / duY^(u)exp 



^Yvm>{u)Pl'm>{r) 



I'm' 



(34) 



with 



Pim(r) = 7jT^2 dr'r'p l 
I'm' 



'm'{r')wimi'm'{r, r') 



(35) 



and 



Pim{r) 



— Y]wm'{r) / dLoYf m {uj)Y Vm ,(Lo)po{r,uj). 
Pb , J 

I'm' 



(36) 



These equations are valid both for the outside and the inside of the sphere or cylinder if 
the r' integrations are taken over the interval [R, oo) or [0, R], respectively. But we note 
that the functions po have completely different forms in these two cases (see Appendix 0). 
It is assumed that the fluid inside a spherical or cylindrical cavity is in equilibrium with a 
particle reservoir at the chemical potential p corresponding to the bulk density p&, which is 
kept fixed when R is varied. For small radii the actual density at the center of the cavity 
may differ from although this effect is certainly numerically neglegible in the examined 
range of radii \R\ > 3. 
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In practice we have truncated all I sums at l max = 10 (l max = 8) for spheres (cylinders) 
and the radial integrals were cut off at a distance C = 2L (£ = 1.5L) from the wall. 
Beyond this distance the profile was assumed to take on its bulk value and corresponding 
asymptotic corrections were added to pi (pi m ) in the vicinity of the cutoff. A step size of 
Ar = 0.02L (Ar = 0.03125L) was used for all functions of r. First the values of iyj 1 i 2 (n, r 2 ) 
(ri,r 2 )) were calculated and stored for all necessary values of n and r 2 and of 
the indices. This step required by far the largest fraction of the computer time. Thereafter 
for a series of bulk densities pb the coefficients $ (vi m ) were determined by a simple Picard 
iteration scheme with retardation. 

In the following a negative (positive) radius R signifies that the wall curves towards 
(away from) the fluid, and z is the distance from the surface. As reduced density we employ 
p* = pDL 2 ; in these units the isotropic-nematic transition takes place at pi ~ 4.2 [|TJ, 
which provides an upper limit for the present approach because in the nematic phase the 
orientational structure does not exhibit the symmetries assumed here. A typical density 
profile p(z, cosO) outside of a sphere is shown in Fig. 0. For z < L/2 orientations with large 
cos# are forbidden so that the profile has a discontinuity along the line cos 9 = x max (z) 
determined in Appendix [B]. When the rods do not interact among each other, i.e., for pi, — > 0, 
all allowed orientations have the same probability. The presence of the steric interaction 
induces a strong increase of the density close to the surface, while there is only a weak 
dependence on cos 9 within the allowed region. Orientations near the discontinuity, where 
one end of the rod touches the wall, are slightly favored. Note that no packing effects are 
visible. These will occur on the much smaller length scale D and presumably only at much 
higher densities where the packing fraction $ ~ pD 2 L is of order unity. The profiles for 
other radii, even for the opposite sign of the curvature, look essentially the same. In the 
latter case there is a very small region close to the surface that is not accessible to any rod 
center. For a cylindrical wall the profiles also depend on the azimuthal angle, but except 
very close to surfaces with negative curvatures this dependence is very weak and a plot of 
p (z, cos#, 0) for any fixed <fi looks very similar to Fig. 0. 

The normalized orientationally averaged number density is defined by 



This function increases for small z up to z = L/2 where it exhibits a cusp and then rapidly 
decreases to its bulk limit 1 which is essentially reached already at z = L. As shown in Fig. |3| 
within the examined range of curvatures (\R\/L > 3) it depends only slightly on R. If 1/R 
is decreased n(z) becomes smaller for z/L < 0.27 and larger for z/L > 0.27. The results for 
the cylinder lie between those for the planar wall (R = oo) and for a sphere with the same 
radius. Due to the finite step size Ar and the steepness of p(z, to) the raw data for n(z) 
exhibit visible kinks at z — nAr for small integers n. These have been removed from Fig. |3| 
by fitting of an appropriate smooth function to the data. 

We define position dependent orientational order parameters as 



For the sphere due to the azimuthal symmetry one has Qi m = for m^fl. The lowest non- 
trivial order parameter Q2o{z) is plotted in Figs. § and |5|. As p(z, u) j (pbTi(z)) ~ 1 for z — » 




(37) 




(38) 
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and in this limit only = tt/2 is allowed it follows that Q2o(z — > 0) = — a/5/87t = —0.3154. 
Negative values of Q 2 o indicate that the rods are preferentially aligned parallel to the surface 
which is of course enforced by the wall. The inter-particle interactions tend to align the rods 
also for z > L/2, where they cannot directly touch the wall, and increase the alignment for 
z < L/2 (see Fig. |J). The alignment is stronger for positive than for negative curvature (see 
Fig. D- 

The biaxiality of the orientational distribution at a cylinder is measured by Q22{z)- 
Positive (negative) values correspond to a preferential orientation perpendicular (parallel) 
to the cylinder axis. If the orientational distribution is sharply peaked at 9 = tt/2 and 
= (or (f) = ix i '2) Q22 takes on its maximum (minimum) value Q22 = ±a/15/327t = 
±0.3863. The results for cylinders of radius R/L = ±5 are displayed in Fig. As expected, 
particles inside a cylinder orient themselves mainly parallel to the cylinder axis, the more the 
stronger the interactions are. With increasing bulk density the decay of Q22 towards the bulk 
value becomes significantly slower, which probably signals the onset of the formation of a 
nematic wetting layer upon approaching the isotropic-nematic transition. This interpetration 
is supported by the fact that the iterations did not converge for pi > 3.5. For this density 
range spontaneous biaxial orientational order has also been predicted at a planar surface 
||. For positive curvature Q22 has the opposite sign and a much lower absolute value. 
The limit z — > for Q22 cannot be determined rigorously in this case because the allowed 
region in ui space does not reduce to a single point so that — > 0) still depends on 
an unknown function of (f). Only in the ideal case p^ — > this function is constant so that 
Q22( z — * 0) = a/5/967t = 0.1288. Here, the effect of the interactions is to increase the 
probability for orientations parallel to the axis so that Q22 may even become negative. 



IV. SURFACE TENSION 

The curvature dependent surface tension "f(R) is that contribution to the grand canonical 
potential which scales with the surface area of the confining wall. In order to determine this 
quantity for positive curvature one has to consider systems of finite size £ in the radial 
direction. However, such systems necessarily contain a second, isotropic liquid - vacuum 
interface generated by the cutoff yielding the corresponding artificial surface contribution 
7 mc . Hence for a sphere we have 



4tt j RV( j R) = lim 



Att 

Q(R, £) -—((R+ £) 3 - i? 3 M - ^(R + £fiL{R + C) 



(39) 



where Uf, = —p (see Eq. (pCf)) is the bulk grand canonical potential density and p the bulk 
pressure. 

The vacuum surface tension can be obtained separately by considering the one-surface 
problem of a sphere of radius R + £ filled completely with an isotropic fluid and in contact 
with the vacuum. In this case the grand canonical potential is the sum of a bulk term and 
a surface term proportional to 7 mc which leads to 

R+C 00 

4n(R + £) 2 {3^ ac (R + £) = ~p 2 b J drr J dr'r'w 00 (r } r'). (40) 

R+C 
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In the deriving this expression we have used the relation 



47rr 



dr'r'w 00 (r,r') = 4t> , 



(41) 



which can be proven using Eqs. fllBD and After some algebra one finally obtains 

oo ( 1 



OO oo 

+ ^Pb J dr i r i -47rrif + y dr 2 r 2 ^2a h 

R L i? /lli2 



Ji)a/ a (r2)tWi lIa (ri,r a ) 



(42) 



Inside a sphere the fluid volume is finite which does not allow to carry out the thermo- 
dynamic limit. Instead we define j s as 



4irR 2 j s (R) = n(R) - —\R\ 3 uj b . 



(43) 

The resulting expression for j(R) is identical to Eq. fl4*2|) , but with all radial integrations 
spanning the interval from to \R\ instead of from R to oo. If one uses the fact that the 
equilibrium profile solves the Euler-Lagrange equation Eq. ( p6|) and that the bulk density 
satisfies the equation lnp&A 3 = (3fi — pbV these results can be cast into the following simpler 
form: 

4irR 2 f3j s (R) = An J drr 2 [l - 2a {r)\ (44) 

~\pI J dr i r ± ~^nvo + J dr2r2^2ai 1 {rx)ai 2 (r 2 )w h i 2 (ri,r 2 ) 

with different integration limits for the outside and inside, as stated above. However, one 
should keep in mind that Eq. ( f4"2"D represents the surface contribution to the density func- 
tional and is minimized by the equilibrium profile whereas Eq. ([H]) only applies to the 
equilibrium solution. In practice the calculation of the surface tension via both formulas 
provides a helpful check of the numerical procedure. 
The corresponding expression for the cylinder is 

2nR 1 c (R) = lim [n(R, C)/H - ir((R + C) 2 - R 2 )u b - 2tt(R + C)^ ac {R + £)] . (45) 

C— >oo 

No truncation in the axial direction is needed because Q is evidently proportional to the 
macroscopic height H . Here we could confirm the analogue of Eq. (|4~l|), i.e., 



dr'r'w QO oo(r,r') = 8tt 2 v , 



(46) 
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only numerically. The resulting expressions for the surface tension, which with the modifi- 
cations quoted above are also valid for the inside of a cylinder, are 

oo 

2nR(3^ c = 2tt J drr^J dup(r, uj) [In 4vrp(r, uj) - 1 -/?//] - p 6 (ln p b \ 3 - 1 - /3/i) | (47) 

oo r oo 
+ ^Pb J dr i r i -2vrf + J dr 2 r 2 ^ Phm 1 (r 1 )p hm2 (r 2 )w h (ri,r 2 



and, at equilibrium, 

oo 

2irR(3j c = 2ti J drrp b 

Jrl 



R 



1 - V A-K/iaoir) 



(48) 



-\ P l / ^ 



y dr x n -27iv + J dr 2 r 2 ^ Pi imi (ri)pi 2m2 (r2)w h (ri,r 2 ) 



We remark that the surface tension depends on the assumed position of the actual surface 
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i.e., on the definition of what is denoted as the volume of the sphere or cylinder, which 
is not uniquely determined. A different choice for this position would alter the volumes and 
surface areas occuring in Eq. (|39| ) and thereby in general lead to a different value of 7. For 
the planar case this change is A7 = pAx, where Ax is the shift in the surface position, while 
more complicated relations emerge for curved surfaces, which may significantly change the 
curvature dependence of what is denoted as the surface tension. On the other hand, experi- 
mentally observable quantities do not depend on this arbitrariness of assigning a particular 
value to the volume of the cavity. For the thin rods we employed the natural definition that 
the defining surface is given by the position of the rod ends at closest approach. But already 
for hard spheres, or for rods of finite thickness, there are at least two possible "natural" 
definitions (see, e.g., Figs. 1 and 14 in Ref. jl8|). 

The surface tension 70 in the ideal gas limit is obtained from Eqs. 



H|) and © by 

neglecting the interaction contributions and by inserting the ideal profiles po(r,uj) from 
Appendix O. This results in 



Po R 2 



drr 2 (l - 2a (r)) 



1 



drr (I — x r 



(r)) 



(49) 



for a sphere (with integration limits for R as described above) and 



Po 



c 



±fdri 

R 
\R\ 

h I drr 



I / dx(p c (r,x) 



l-fdx(l-^ c (r,x)) 



R>0 



R<0 



(50) 



for a cylinder. This limit has been discussed extensively by Yaman et al. who found 

the surprising result fi^o/po = L/A for rods outside any convex body, while this expression 
is modified inside a sphere or a cylinder, so that 
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^ (51) 



Po If, R > o 

and 

po "If, R>o ■ {b2) 

We have confirmed these results (analytically for the sphere, numerically for the cylinder) 
by using Eqs. ([49]) and (|50|) . The most interesting aspect of these findings is that the surface 
tension is not analytical at 1/R = which may lead to unexpected behavior of membranes 
immersed in colloidal rod suspensions. 

We have determined the surface tension in the presence of the interparticle interactions 
for a series of bulk densities and radii. The results are shown in Fig. [7] as function of p& for 
fixed R and in Fig. |8] as function of R for fixed p&. For each data point four (spheres) or three 
(cylinders) numerical calculations with different values of the cutoff l max were performed. 
The results were extrapolated to l max = oo using a fit function linear or quadratic in l/l max . 
The differences betweeen the values at the largest l max and the extrapolation become consid- 
erable [A(flj) / (pbL) ~ 0.02] at large densities. From comparison of the results obtained by 
quadratic and linear extrapolation we estimate the error of (3 r y/{pbL) to be 0.01 for pi = 4 
but only 0.001 for pi < 2. Finally we have interpolated smoothly between the 13 data points 
taken for each radius. In the case of a cylinder at the highest density pi = 3 the result 
obtained separately for the special case of a planar wall lies slightly (by 0.003 in the units 
used here) above the almost linear curve through the other points. A possible explanation 
for this observation is that at the planar wall a uniaxial orientational distribution has been 
assumed while the actual equilibrium profile might exhibit a small spontaneous biaxiality 
as found at the cylindrical walls. Therefore the planar wall result has not been used for the 
interpolation scheme in this case. 

Figure [I] shows the surface tension divided by the density to allow for a better comparison 
with the ideal rods results that predict a density independent constant for this ratio. In all 
cases the interaction significantly increases this quantity, by up to 50% for the intermediate 
density pi = 2. In the upper density range saturation or the formation of a maximum are 
observed. The results for a planar wall are in good agreement with those obtained by Mao et 



al. [O]], who effectively used the same theory but a different numerical method. For almost 
all densities the surface tension is higher for negative than for positive curvature, in contrast 
to the behavior at pb — > 0. In the latter case 7 is constant for R > 0; the slight R dependence 
for R < is hardly visible on the scale of Fig. ||. On the other hand, for densities of the 
order of 1 in reduced units the dependence on R is dominated by a term linear in 1/R that 
is absent in the ideal limit. The dependence on R becomes stronger and less linear with 
increasing density, especially for the spherical case. 

Due to the limited number of radii for which calculations were performed we are not able 
to decide whether the small discontinuity of the second derivative of j(l/R) at 1/R = 
that occurs in the ideal limit persists also at finite bulk densities. However, our data do not 
preclude this possibility. It is commonly assumed that the surface free energy density for a 
general surface with principal curvatures 1/R\ and I/R2 has the Helfrich form |TP| 

7 (l/i?i, I/R2) =1 P + 2k (hi- + h) _ Co ) + =1_ + . . . (53) 



2 Ri R2 J RiR 
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with the stiffness coefficients k and R and the spontaneous curvature cq. Inter alia, this 
form predicts that the contribution to j(l/R) linear in 1/R for a sphere is twice that for a 
cylinder. From our numerical results we found that this relation is approximately fulfilled at 
low densities, but there are substantial deviations from it at higher densities. Moreover, the 
Helfrich expression cannot be applicable for all signs of the curvatures already in the ideal 
limit due to the aforementioned non-analyticity of 7 S (l/i?) and ^(l/R). 

We conclude that for phenomena for which the curvature dependence of the surface ten- 
sion matters the non-interacting approximation is qualitatively wrong already at relatively 
low densities (compared to the transition density to the nematic phase) and at any fixed 
curvature it produces rather large quantitative errors in the absolute value of 7. 

V. SUMMARY 

For an isotropic fluid of needle-like hard sphero cylinders of length L near hard spherical 
or cylindrical walls we have obtained the following main results: 

1. Due to the interparticle interactions the probability of finding a particle with a given 
orientation is strongly increased close to the wall when compared to the bulk fluid. At 
a given distance z of the center of mass from the wall orientations for which one end 
of the rod touches the wall are most favorable (see Fig. |2[). 

2. Since the range of accessible orientations decreases when the particle approaches the 
wall the orientationally averaged density vanishes for z — > 0. It exhibits a cusp at 
z = L/2 where the rods lose contact with the surface (Fig. H). 

3. The parallel alignment favored by the surface decays more slowly when the bulk density 
is increased (Fig. |) and is stronger if the wall curves towards the fluid instead of away 
from it (Fig. |5|). 

4. A cylindrical wall curving towards the fluid induces biaxial orientational order with 
preferential alignment parallel to the cylinder axis. With increasing bulk density a 
nematic wetting layer develops in this case (see Fig. |6|). 

5. The density and curvature dependences of the wall-fluid surface tension are shown in 
Figs. |7] and In contrast to the results for non-interacting particles given by Eqs. Q5TD 
and (^) the surface tension decreases with increasing curvature and exhibits a linear 
behavior around the planar limit which leads to a spontaneous curvature of a membrane 
away from a fluid of rod-like (colloidal) particles. 



APPENDIX A: DETERMINATION OF THE EXPANSION COEFFICIENTS OF 

THE MAYER FUNCTION 

Equation (0) serves as the starting point which defines the expansion coefficients fui 2m (r). 
As far as the azimuthal angles are concerned the integrand depends only on 12 = <p2 — 4>i 
because 
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Y lm {9,<j>) = k lm P lm (cos6)e- m t (Al) 

where the coefficients ki m relating the spherical harmonics to the associated Legendre func- 
tions Pi m are given by 

km = (-1)™ ( 2 4^ (1 + ^)1 ) ' fOT m ^ °' = (- 1 ) m ^- ( A2 ) 

Therefore after the substitutions 4> s = \{4>i + 4>2) and i2 = 4>2 — 4>i the integral over S 
renders a factor 2ir. In Eq. ([?[) we now perform the integrations in the following order: 

dcos^i / dcos6> 2 / d4) l2 f(rj l7 6 2 , 0i2)-P;im(cos^i) (A3) 

1 J-l Jo 

xP hm (cos6 2 )e im ^. 

In the limit D/L — > the two rods overlap only in a small range A0 12 around 0i2 = 
or 12 = 7T, approximately given by A0 12 = 2D/p where p is the distance between the 
intersection point and the line joining the centers of the rods, as shown in Fig. |9|. From this 
figure one easily derives p = r/(cot^i — cot #2); based on appropiately modified figures one 
finds that this expression is also valid if Q\ > 71/ 2 or 62 < it/ 2. Therefore it is sufficient 
to replace the (fru integration by the factor A<pi2, to replace the integrand by its value at 
012 = or 7T, and to examine only the overlap of two infinitely thin rods in the xz plane. In 
integrations over the full spatial angles Cji we assign an intrinsic directionality to the rods 
such that their "front ends" point into the direction (^j,0i) and their "rear end" into the 
direction (it — 0i,2n — 0j). In order to calculate fi^m it is sufficient to integrate over the 
configurations with overlap of the two front ends, where necessarily 12 — 0. One can easily 
show by appropriate substitutions that the other three possibilities for overlap yield exactly 
the same contribution to the total integral, so that one finds 



D f 1 

fhhm\r) = -16irk hm ki 2 m— / dx 1 Pi im (x 1 ) (A4) 



D f 1 

x max{ r / L,X±) 



x J dx 2 Pi 2m (x2) 



£l X 2 



.1 ■-) 



Here x min and x max denote the smallest and largest value of x 2 = cos 9 2 for given X\ = cos Q\ 
and r/L for which the front halves of the rods overlap, as shown for an example in the 
inset of Fig. They can be determined by tedious but straightforward geometry, which 
yields two basic formulas for the x value at touching, depending on whether the end of rod 
1 touches rod 2 or vice versa: 



xf ] (r', x) = — xf' ±} (r', x) = 2r'(l - x 2 ) ± xy/l - 4r' 2 (l - x 2 ) (A5) 



yl — 2r'x + r 

with r' = r/L. As illustrated in Fig. [T^ the following regions in the (r',x) plane must be 
distinguished: 
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A, for < r' < 1/2 and r' < x < 1: x m j n = and x mas = 1, 

B, for 1/2 < r' < 1 and r' < x < 1: x min = and x maa; = x^ 2,+ ' ) , 

C, for 1/2 < r' < l/y/2 and ry/l- l/4r' 2 < a; < r': x min = xf~ ] and x max = xf' +) , 

D, for < r' < 1/2 and — 1 < x < r'\ x m i n = x^'~^ and x max = 1. 



For given I and m the integration over x 2 in Eq. ( |A4| ) can be carried out analytically. 
In the remaining numerical integration special care must be taken for r' < 1/2 due to the 
square root singularity at X\ — 1. All coefficients fi x i 2m are of the order of D/L and vanish 
for r/L > 1. They diverge for r/L — > 0, but they appear only in the product r/ iliam (r) 
(see Eqs. (pT3|), ([16]), and (|19D) which is finite in this limit. For two finite values of D/L the 
coefficients fi^m have been calculated by Moore and McMullen 0. 

A useful check of the numerical results is obtained from the observation that the excluded 
volume t> e;z . (cos 7) for fixed angle 7 between the particle axes is related to the Mayer function 
via v ex (cos'y) = — f d 3 r 12 /(ri 2 , u>i, u 2 ) from which one derives 

dcos 7 PKcos7)^(cos7) = — I ^(-l) m+1 / drr 2 f Um (r). (A6) 

The left hand side can easily be determined using the well-known result (see, e.g., Ref. [0) 
fez (cos 7) = 2DL 2 \ sin 7 1 + 0(D 2 L) and thus provides a sum rule for the second moments of 
the expansion coefficients. Our numerical results passed this check. 



APPENDIX B: DENSITY PROFILES FOR NON-INTERACTING RODS 

In the ideal limit p{r,u) adopts the constant value po/4vr for orientations to that are 
allowed by the presence of the hard wall and vanishes otherwise. Thus it is sufficient to 
determine the limiting orientations for which the rod just touches the wall. At a planar 
surface the maximum allowed value for x = cos 8 is x max = 2z/L for z < L/2 while the rod 
cannot touch the wall for z > L/2 so that x max (z > L/2) = 1. At all surfaces the minimum 
value of x clearly of the head-tail symmetry. Thus in the following it is 

sufficient to consider positive x. 



1. Outside a sphere 

As illustrated in Fig. [TT|(a), there are two different ways how a rod can touch the outside 
of a sphere: when the rod is sufficiently far from the surface its end touches the wall upon 
rotation, whereas when it is close to the wall at contact it will touch it tangentially. The 
crossover between these two regions takes place at r c = \J R 2 + L 2 /4 and straightforward 
geometrical reasoning yields 

1, r>R + L/2 

./■„„,, = < (L 2 /4 -R 2 + r 2 )/(rL), r c <r<R + L/2 (Bl) 
\/r 2 — R 2 /r, R < r < r c 
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2. Inside a sphere 



In this case there is a minimum distance r c = \J R 2 — L 2 /A from the surface beyond 
which all orientations are forbidden. In the accessible region \R\ — L/2 < r < r c one finds 
[see Fig. 0(a)] 



R 2 - L 2 /A 



r 



2 



•£"tnax r • V / 

rL 



3. Outside a cylinder 

For a cylinder the profile depends in addition on the azimuthal angle which we always 
measure from the axis that is perpendicular to both the cylinder axis and the surface normal. 
Thus here we have to determine the range 1$ of allowed values of for fixed r and cos (9. 
If the rod touches the cylinder at the angle C G [0, tt/2] then due to symmetry one has 
1^ = [0, C ] U [n — C , 7r + C ] U [2n — C , 2ir]. With the help of Fig. |TT](b) one finds for contact 
between the rod end and the surface 

cos 0c = — — B3 

L sin 6 

whereas for tangential contact one has 

R 2 

cos 0c = cot 9 (B4) 
yr 2 — R 2 

By determining the transitions between these two well as those to the ranges 1^ = 

and 1$ = [0, 2n] we are led to distinguish the following cases (see Fig. |12"D : 



A, R < r < R + L/2 and < cos9 < 2{r - R)/L: cos0 c = 0, 



B, r c = ^JR 2 + L 2 /4 < r < R + L/2 and 2(r -R)/L< cos$ < (L 2 /4 - R 2 + r 2 )/(Lr): 
cos 0c as given by Eq. (|B~3|), 

C, R <r <r c and 2(r — R)/L < cos6 < 2{r 2 — R 2 )/(Lr): cos0 c as given by Eq. (|B3|), 



D, R < r < r c and 2(r 2 — R 2 )/(Lr) < cos 6* < \Jr 2 — R 2 /r: cos0 c as given by Eq. (|B 

In the remaining region within R < r < R + L /2 all values are forbidden, while obviously 
for r > i? + L/2 all orientations are allowed. 



4. Inside a cylinder 

In contrast to the previous case, here the accessible range is centered around the 
cylinder axis, i.e., it has the form 1$ = [0 C , 7r — C ] U [n + 4> c ,2ti — C ]. Since tangential 
contact is not possible, the classification is a little bit simpler. In the region of interest 
\R\ — L/2 < r < \R\ we obtain: 
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A, for cos 9 > 2(r - \R\)/L\ cos0 c = w/2, 



B, for r > r c and cos# < 2(r — \R\)/L or for r < r c and (i? 2 — L 2 /4 — r 2 )/(Lr) < cos^ < 



C, for r < r c and cos6> < (i? 2 — L 2 /4 — r 2 )/(Lr): cos0 c = 0. 

In Fig. [H| we compare the accessible orientational space in these four cases for the same 
radius \R\fL = 3 and the same distance from the surface z/L = 0.2. Naturally this space is 
largest outside a sphere and smallest inside a sphere. One also notices that for a cylinder the 
<j) dependence is actually restricted to a rather small range of values for cos 9 while for most 
polar angles 6 either none or all azimuthal angles are allowed. We emphasize again that the 
forbidden regions are also strictly forbidden for interacting rods whose profile is no longer 
constant within the allowed region. 



2{r-\R\)/L: 



COS (J), 




(B5) 



'c 



L sin 6 
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FIGURES 



FIG. 1. The system under consideration consists of a fluid of monodisperse hard spherocylinders 
of diameter D and length L in contact with a spherical or cylindrical hard wall of radius R. We 
assume that the orientational distribution only depends on the normal distance z and, for a spherical 
wall, the angle 9 of the particle axis with respect to the surface normal. During the evaluation of 
the excess free energy the particle orientations are described in three different reference frames: the 
frame x'y'z' fixed in space, the frame xyz fixed by the local normal direction, and the interparticle 
frame xyz. The rotation between the latter two is described by the Euler angles ipi, rji, and Xi ( see 
Eq. (||)). The different sizes of the spherocylinders indicate that the particles typically do not lie in 
the x'y' plane and thus only their projection onto that plane is shown. For reasons of clarity only 
the polar angles 9', 9, and 9 are shown, but not the corresponding azimuthal angles cj>', 4>, and (f>. 

FIG. 2. Full density profile p*(z, cos 9) = p(z, cos 9)DL 2 outside a sphere of radius R/L = 3 
for the bulk density p\ = pbDL 2 = 2. At small distances z from the wall large values of cos 9 are 
forbidden due to overlap. Therefore the profile is exactly zero behind this "step". In the allowed 
region the most prominent feature induced by the interaction between the rods is the strong increase 
of the density at small distances z. For fixed z orientations close to the step, i.e., those with one 
end of the rod touching the wall, are favored. If the interparticle interactions were neglected the 
profile would be constant in the whole accessible region. 

FIG. 3. Normalized orientationally averaged density n{z) for fluids in contact with planar, 
cylindrical, and spherical walls of positive and negative curvature for a fixed bulk density pi = 2. 
The inset shows the behavior in the vicinity of the cusp which occurs at z = L/2 and is followed by 
a rapid decay towards the bulk limit n{z) = 1. All curves have about the same value at z/L ~ 0.27 
but they do not intersect exactly at one point. For better visibility only two curves are shown in 
the main part of the figure. The remaining profiles lie in between these two. Close to the wall the 
number density is larger for positive curvature. Inside a sphere n(z) = for very small z because 
the centers of the rods cannot come arbitrarily close to the wall. Small kinks in the small z range 
that are caused by the numerical discretization were removed by fitting a smooth curve to the raw 
data. 

FIG. 4. Uniaxial nematic order parameter Q20 as defined in Eq. ( |3"8| ) for a fluid outside a 
sphere of radius R/L = 5 for different densities. The tendency for parallel orientations, as indicated 
by negative values of Q2O1 becomes more pronounced and longer ranged with increasing density. 
Q2o{z = 0) = -0.3154. 

FIG. 5. Order parameter Q20 at fixed density p\ = 2 for different wall geometries. Parallel 
orientations, i.e., negative values of Q20 are more favored by negative curvatures. 

FIG. 6. Biaxial order parameter Q22 (see Eq. (|3S|)) for cylindrical walls of positive and negative 
curvature at different bulk densities. The most probable orientations for positive and negative 
values of Q22 are indicated in the sketches. The strong increase of the decay length with density 
for R/L = —5 signals the formation of a nematic wetting layer. 
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FIG. 7. Density dependence of the surface tension for different wall geometries and curvatures. 
In the ideal limit fi^f/p^L takes on the density independent values 0.25 for 1/R > (cylinders, 
spheres, and planar wall), and 0.2492 (0.2497) for a sphere (cylinder) with R/L = —5. The differ- 
ences between these three values are not visible on the scale of the figure. 

FIG. 8. Curvature dependence of the surface tension for (a) spheres and (b) cylinders at various 
densities. In the limit of non-interacting rods (pb = 0) there is a slight decrease of 7 for negative 
curvature, which is hardly visible on this scale. Taking into account the interparticle interactions 
leads to a substantial enhancement of the curvature dependence already for pf = 1. We have 
interpolated smoothly between the 7 data points calculated for each curve. 

FIG. 9. Geometry of two overlapping rods in the limit D/L — > 0. Only one half of each rod is 
shown. Their centers lie on the z axis and have a distance r. In order to overlap both rods must lie 
approximately in the same plane (a). Part (b) shows a vertical projection from which the allowed 
azimuthal range lS.<p\2 = 2D/p can be determined. Rod 2 is drawn in the two positions for which 
it just touches rod 1. 

FIG. 10. Illustration of the cases that must be distinguished for the determination of the 
minimum (x m i n = cos9^ iax ^ and maximum {x max = cos9 2 ntn ) value of X2 so that two rods overlap 
for given x\ = cos 9\ and r. The appropriate expressions for x m i n and x max for the regions A, B, 
C, and D are given in the main text. In the unlabeled region containing the inset overlap is not 
possible. The inset shows an example from region B where the end of rod 2 touches rod 1 at the 
minimum angle and vice versa at the maximum angle. 



FIG. 11. Geometries for the determination of the allowed orientations at curved walls, (a) Out- 
side a spherical wall the rod touches the wall with its end if r > r c (rod 1), but tangentially if r < r c 
(rod 2). Inside a sphere only end contact can occur (rod 3). From this figure Eqs. (|BlD and (|B2|) can 
be derived, (b) Projection of the corresponding problem for a cylindrical wall onto the plane per- 
pendicular to the cylinder axis. Equations (|B3p and Q can be derived using xp = ^ sin 9 cos cf> c 
and zp = j cos 9. 



FIG. 12. In calculating the accessible orientational space for a rod outside a cylinder different 
expressions, presented in the main text, apply in the regions labeled A, B, C, and D for the 
maximum azimuthal angle <p c at given values of x = cos 9 and r. The figure corresponds to the 
case R/L = 1.3, but its topology is the same for all radii. Only for configurations corresponding to 
region D the rod touches the cylinder tangentially. For larger values of R as actually used in the 
calculations presented above r c is closer to R and the regions C and D are much smaller. 

FIG. 13. Accessible orientational space for a rod near walls of different geometries at a fixed 
distance z/L = 0.2 and radius \R\/L = 3. Spatial angles to the left of the lines are allowed, those to 
the right are forbidden by the presence of the wall. There is no dependence on the azimuthal angle 
(j) for spherical walls. For a cylinder at (j) = ((f) = tt/2) the wall is effectively spherical (planar) 
which explains the common end points of the various lines. 
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